Possibility of superconductivity due to electron-phonon interaction in graphene 
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We discuss the possibility of superconductivity in graphene taking into account both electron- 
phonon and electron-electron Coulomb interactions. The analysis is carried out assuming that the 
Fermi energy is far away from the Dirac points, such that the density of the particles (electrons 
or holes) is high. We derive proper Eliashberg equations that allow us to estimate the critical 
superconducting temperature. The most favorable is pairing of electrons belonging to different 
valleys in the spectrum. By using values of electron-phonon coupling estimated in other publications 
we obtain the critical temperature Tc as a function of the electron (hole) density. This temperature 
can reach the order of 10 K at the Fermi energy of order 1 — 2 eV. We show that the dependence of 
the intervalley pairing on the impurity concentration should be weak. 

PACS numbers: 74.70.Wz , 74.62. En, 74.25.Bt 



I. INTRODUCTION 

Since its experimental discovery in 2004^ and first ob- 
servations of unusual properties, ■^'^ graphene has gained 
a lot of experimental and theoretical attraction. In the 
last decade, thousands of articles devoted to the study 
of graphene have appeared. Possessing novel electro- 
mechanical properties, graphene is a promising material 
for electronic devices. The linear Dirac- type electron 
spectrum makes graphene very interesting from the the- 
oretical point of view and many interesting effects have 
been predicted and observedi^ 

However, still there are open questions on fundamental 
electronic properties of graphene, and one of the most in- 
teresting ones concerns a possibility of superconductivity. 
Graphene is a good conductor unless the Fermi energy is 
too close to the Dirac points and the electron-phonon 
coupling in graphene is not very weak. Therefore, al- 
though the superconductivity in graphene has not been 
observed, it is not clear why one should discard the pos- 
sibility of this phenomenon. 

Superconductivity can be induced in graphene by su- 
perconducting contacts due to the proximity effect;^ but 
can it be obtained in a "natural way"? What should one 
do in order to obtain a considerable value of the super- 
conducting critical temperature Tc? What type of the 
superconductivity and what structure of the order pa- 
rameter could one expect? 

In the last years, there have been various attempts to 
answer these questions. Due to the special type of the 
spectrum, the main attention has been devoted to inves- 
tigating the possibility of unusual types of the supercon- 
ducting pairing. Superconducting pairing mediated by 
conventional electron-phonon or electron-plasmon cou- 
pling was considered in Ref. [6| with a conclusion that, 
in addition to the conventional s-wave pairing, an exotic 
p + ip state is possible. Superconducting properties of 
Dirac electrons in graphene were considered in Ref. 
within the conventional BCS approach. In these pub- 
lications, the main emphasis was put on the study of 
properties of unusual superconducting pairing for a small 



electron density. 

It is clear that one can expect very interesting new 
properties of the superconductivity in the vicinity of the 
Dirac points. However, this region is least favorable for 
the existence of superconductivity due to the very low 
density of states, and one should tune the Fermi energy 
away from the Dirac points in order to have a hope to 
obtain superconductivity. 

By doping graphene by various combinations of K and 
Ca, the authors of Ref. [sl were able to shift the Fermi 
energy far away the Dirac points and even to put it in the 
vicinity of the van Hove singularity (VHS). Another ex- 
perimental method based on the use of electrolytic gates^ 
allowed the authors to tune continuously the electron 
density up to values n — 4 x 10^*cm^^, which is ap- 
parently not very far away from the VHS. These experi- 
mental works have demonstrated that one can achieve a 
ultrahigh electron density and this makes observation of 
superconductivity in graphene considerably more realis- 
tic. At the same time, transport measurements were not 
carried out in Ref. fs'l and the superconductivity has not 
been seen in Ref. [9] for temperatures T > 1 K. 

Although the superconductivity has not been observed 
yet, theoretical considerations^''"^^ predict superconduc- 
tivity at the VHS even for a repulsive electron-electron 
interaction. In this case, the superconductivity is ex- 
pected to have an unconventional symmetry of the order 
parameter. No doubt, an experimental observation of 
the superconductivity at the VHS would be of a great 
interest, but one should be able to tune the Fermi en- 
ergy exactly to the singularity. Disorder may also play a 
destructive role in formation of such a superconductivity. 

Therefore, it would still be important to investigate 
theoretically the possibility of a superconducting pair- 
ing due to the conventional electron-phonon pairing far 
away from the Dirac point, but, at the same time, not 
in the vicinity of the VHS. Such a study implies using 
conventional schemes of computing the superconducting 
transition temperature. Then, one should check the sta- 
bility of the pairing against disorder in the system, clarify 
the dependence of the transition temperature Tc, etc. 
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In several publications, models with an electron- 
phonon interaction have been considered. The authors 
of Ref. discussed the valley structure of the or- 
der parameter using an electron-phonon model with- 
out the electron-electron Coulomb interaction. They ar- 
gued that there might be a superconducting instabil- 
ity in highly doped graphene, while the valley structure 
depends on the parameters of the electron-phonon in- 
teraction. Superconductivity in hydrogenized graphene 
(graphane) has been considered as wellji^ In this sys- 
tem, a model based on electron-phonon interaction was 
used and a transition temperature of 90 K in p-doped 
graphane was predicted. At the same time, it is clear 
that taking into account the Coulomb interaction is very 
important because it can in principle be even stronger 
than the electron-phonon interaction. Moreover, consid- 
ering the latter in the weak coupling limit can also lead 
to incorrect predictions. 

In this paper, we use a generic model including both 
the electron-phonon and electron-electron Coulomb in- 
teractions. Using Eliashberg-type equations,^—"— we de- 
rive an expression for the transition temperature de- 
termined by the electron and phonon interactions in 
graphene. By using experimental values and results of 
numerical calculations for the interactions obtained for 
the normal state, we estimate the transition temperature 
Tc and conclude that the superconductivity is possible 
with Tc of the order of several Kelvin. Most favorable is 
a singlet pairing between different valleys. We show that 
such a pairing is not very sensitive to disorder. 

The paper is organized as follows. In Sec. II, we formu- 
late the Hamiltonian of quasiparticles in doped graphene 
describing the interaction of quasiparticles with phonons 
and the Coulomb interaction. In Sec. Ill, we consider the 
electron pairing by deriving and solving the Eliashberg 
equations. The effects of impurity scattering are consid- 
ered in Sec. IV. In Sec. V, estimates are made, and the 
concluding Sec. VI is devoted to discussions. 



II. MODEL HAMILTONIAN 

In this section, we introduce the Hamiltonian of inter- 
acting quasiparticles in a single layer of graphene. For 
the undoped system, the Fermi surface consists of two 
nonequivalent points K and K' = — K, called the Dirac 
points.** The quasiparticles around these points have a 
linear spectrum. This approximation remains valid up 
to quite high energies. We consider in this work doping 
levels corresponding to the Fermi energy ep of the order 
of leV. In other words, we consider the case when the 
Fermi energy is sufficiently far away from both the Dirac 
point and the VHS. In this region of parameters the spec- 
trum consists of two well resolved valleys and each valley 
is a cone. As the Fermi energy is far away from the Dirac 
points one does not need to account for effects specific for 
the Dirac equation. 

The two valleys are numerated by a variable called 



isospin. Due to the electron-phonon and electron- 
electron interactions there are matrix elements of the 
Hamiltonian mixing these two valleys. As we do not in- 
vestigate properties of the system near the Dirac point 
we do not use the Dirac-type representation of the 
Schrodinger equation. 

The Hamiltonian describing the electron-phonon sys- 
tem can be written in a general form 

H = Hq + He.ph + He^e, (2.1) 

where 

Ho = J2ie{p)-fi)cl^Cp^, (2.2) 

p,(T 

is the operator of the kinetic energy, e (p) is the spec- 
trum of the non-interacting electrons, /i is the chemical 
potential (Fermi energy at low temperatures), a is the 
spin index, Cp^a {cp is the electron annihilation (cre- 
ation) operator for the electron with momentum p and 
spin (T, He^ph stands for the electron-phonon interaction, 
and He,e-ioT the electron-electron one. 

A. Electron-phonon interaction 

The Hamiltonian He^ph describing the interaction be- 
tween the electrons and phonons can be written in a gen- 
eral form asiSii^ 

He,ph = 9p,ci.J *qjcJ,+q,^Cp,<„ (2.3) 

where 5p,qj is the electron-phonon coupling function and 
^qj- is the phonon field operator for the mode j. As 
usual, the most important contributions to the thermo- 
dynamics are expected from the vicinity of the Fermi 
surface that consists in the case involved of two circles. 

In order to distinguish explicitly between the valleys 
we write ck+p,ct = oi.p.o- and CK'-i-p,cr = a2,p,tT for quasi- 
particles in the vicinity of the Fermi circles of the two 
valleys 1 and 2, where Oc.p.o-, a = 1, 2 are fermionic anni- 
hilation operators with the momentum p measured from 
the Dirac point of the a- valley. 

Using these notations we write the electron-phonon in- 
teraction He.ph, Eq. (|2.3p . in the form 

2 

He,pH = E E C *q>l.P+q,.«^^P.- (2-4) 
a, 0=1 p,q,J,o- 

In Eq. (|2.4p the coupling constants g"^ and the phonon 
field operators are related to those in Eq. (|2.3p as 

5'p!qj =.9p-K,qj, 5p^qj =5p-K',q,j, (2.5) 

12 _ 21 _ 

5p,qj — .9p-K',q-Q,j, .9p,qj — ffp-K.q+QJ, 

$11. — <J> <I>22 ^ 

qj ^qj' ^qj ^qj' 

*qj = *q-Q,J, K]j = *q+Qj' 
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where Q = K' — K is the vector connecting the Dirac 
points. The phonon fields <i> are real in the coordinate 
representation and therefore we obtain the following re- 
lations 



21 



-qj V^qj 



$22 . ^ ('$22 

-qj \ 1,3/ 

$21 . ^ f$12.)* 

-qj V^qjy 



(2.6) 



In Eq. p.4p the sum is taken over such p and q that 
both p and p + q are in the vicinity of the Fermi cir- 
cle. The terms with a = /? in the Hamiltonian H^^ph, 
Eq. (|2.4p describes the intravalley scattering of electrons 
by the phonons, while the terms with a ^ /3 stand for 
the intervalley scattering . Formally, we can speak of an 
isospin dependence of the electron-phonon interaction. 

The bare Hamiltonian Hq. Eq. (|2.2p . takes in these 
notations the form 



(2.7) 



(2., 



Q— 1 p,cr 

where the energy equals 

Cp = I'o IpI - 

and vq « 10® cm s~^ is the Fermi velocity. 

Equations (|2.4p - (|2.7p specify how the presence of the 
two valleys can be written in terms of the isospin. The 
coupling constants g°^^ are different for the intravalley 
and intervalley scattering (equal or unequal a, /3, respec- 
tively). In graphene, two two-phonon peaks are seen 
in the Raman spectrum (see, e.g., Refs. [2O and [2lh : 
the D* peak near 2a;Ai = 2650 cm~^ corresponding to 
the scalar Ai optical phonons and the G* peak near 
2ijjE2 = 3250 cm"-'^ corresponding to twofold-degenerate 
pseudovector E2 phonon mode. The E2 mode is respon- 
sible for intravalley scattering, while the scalar Ai optical 
mode leads to the intervalley scattering^^^ 



B. Coulomb interaction 

The operator of the electron-electron interaction Hg.e 
in Eq. (j2.ip can be written in the standard form 



He.e — 



2 ^ 
pp'q 



S+qS' 



(2.9) 



where is the matrix element of the Coulomb potential 
in two dimensions. We assume here that corrections to 
the bare Coulomb potential have already been calculated 
and use therefore for the static screened Coulomb in- 
teraction. Due to a specific form of the wave functions 
of graphene leading to a suppression of the backscatter- 
ing, the function differs from the conventional Fourier 
transform of the Coulomb interaction and can be written 
in the form 



2^e2 



K|q|e(q) 



7(q) 



(2.10) 



where k is the dielectric constant of the substrate, e(q) = 
1 -I- ^y^n (q) is the static dielectric permeability of the 
electron gas and 



7(q) = \ (1 +C0S(/)q) 



(2.11) 



with 0q = arctan ^ is the form factor accounting for 
the absence of baclcscattering in graphene. The static 
polarizability 11 (q) is given by23. 



n(q) 



for q < 2pp 

^i^"' (^) + for , > 2p, 

(2.12) 



where 



(2.13) 



is the density of states per one spin direction and one 
valley and is the Fermi momentum. The factor 4 in 
Eq. (|2.12p accounts for the number of the spin and isospin 
directions. 

Again, we can use the isospin formulation to bring the 
operator H^.e to a more convenient form 

^-.■^^JE E ^q"''4,p+q4,p'-q«".P'«,9,P. (2-14) 



where 



a,l3 P,P':q 



yll y22 Y 



(2.15) 



■q 'q-W 'q 

and the sum in Eq. (|2.14p includes such momenta p, p', 
and q that both p + q and p' — q are in the vicinity of 
the Fermi circle. 

In the next section we will derive Eliashberg equations 
for the model described by Eqs. I^Ul . (|^ . (j^ . and 
(EH. 



III. ELIASHBERG EQUATIONS 

In order to describe the electron pairing mediated 
by electron-phonon interaction, we derive the Eliashberg 
equation o^^^^'^^'^'^^'^^ for the system under consideration 
using conventional methods of quantum field theory.— 
We introduce the imaginary time normal and anomalous 
Green functions that are 4x4 matrices in the valley and 
spin space 

Gp(r) ^ -(T,^p(0)^p(r)) F^ir) ^ (TW3-p(0)V'^(r)), 

(3.1) 

where the field operators tjjp (r) are 4-components vec- 
tors having as components the operators aa.p.a in the 
Heisenberg representation. 
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In the isospin representation the gap function Ap^^^ 
entering the Gorkov equation is a 4 x 4 matrix. 

Using the normal Gp^e„ ^'Hd anomalous -Fp Green 
functions, we consider the Dyson equations containing 
both normal and anomalous self-energy parts. Using this 
matrix representation, we derive the Gorkov equations in 
the generalized form 

(«^„£„ - Cp) Gp,e„ + Z„Ap,,„F+,,^ = 1 (3.2) 
(zZ„£„ + ^p) F+^^ + Z„A+,^Gp,e„ = 0, 

where £„ = ttT (2n + 1) is the fermionic Matsubara fre- 
quency, and the symbol " -t- " stands for the Hermitian 
conjugation of the 4x4 matrices. 

According to the Eliashberg^'^ theory developed for 
an arbitrary value of the coupling in the limit ujjj <^ 
/i, where loo is the Debye frequency, one has to take 
into account normal contributions of the self-energy to 
the Green functions, but neglect the renormalization 
of the vertices. As we assume that both the electron- 
phonon and electron-electron interactions are not neces- 
sarily weak, one should introduce the factor Z„ coming 
from the normal self-energy. This factor renormalizes 
the coefficient in front of the frequency, but the renor- 
malization of the coefficient for ^p is neglected. This 
non-equivalence originates from the fact that the renor- 
malization of the coefficient in front of the frequency is 
proportional to w^^, while corrections to the coefficient 
in front of ^p are proportional to The factor Z„ in 
front of Ap is written for a convenience. 



By solving Eqs. (|3.2p . we write the anomalous Green 
function F"*" as 



Z,2,A+,„Ap,„ 



(3.3) 

Explicit calculations show that, in the absence of mag- 
netic interactions in the system, the most favorable is 
singlet pairing. Therefore, we do not present here calcu- 
lations for the general case and consider only the singlet 
pairing. At the same time, the structure of the gap func- 
tion Ap^e^ can be non-trivial due to the presence of two 
valleys. In our representation using the isospin, Ap^^^ 
is a 4 X 4 matrix in both spin and isospin space and we 
write it as 



A. 



P,£r 



M" A12 

\^p,e„ ^P,e„ 



(3.4) 



where (T2 is the second Pauli matrix. The fact that only 
this matrix enters the gap function in the spin space is 
standard for the singlet pairing. Exchanging the elec- 
trons the gap function must change the sign. As the 
matrix a2 is antisymmetric, one comes to the relation 



'-'p,e„ ^ — p,e„ 



(3.5) 



Considering the triplet order parameter, one would 
obtain instead of the symmetric relation for Ap^^ 
[Eq. p.Sp ] an antisymmetric one. 

As the spectra of the valleys are identical, we can con- 
sider a simpler form of Eq. (|3.4p taking 



— Ao,p,e„ — ^o,-p,e„, 



= A^ 



= AQ,p,e„. 



(3.6) 



The gap Ao,p,e„ corresponds to the intravalley pairing 
and Aq.p.£^ to the intervalley one. 

By substituting Eqs. p. 41) and p.6p into p.3p . we write 
the solution for the anomalous Green function F, 



P,£r 



as 



Z„(Pp+,„+Pp-;,J 
2 



Ao,p,£„ 




Aq-p^e,. 


A0,p,£, 


Aq -p,e„ 


Ao,p,£, 


Ao,p,e„ 


AQ,p,e 



5 10-2 

(3.7) 



where the function Pps„ equals 



[Z,,en? +e, + Zl |Ao,p,e„ ± Aq,p,, 



(3.8) 



Equation (|3.7p should be complemented by a self- 
consistency equation, which is actually a matrix equation 
in the isospin space. Writing separately the anomalous 
and normal parts of the self-energy, we come to the fol- 
lowing equations: 
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where Dj is the phonon Green function for the polar- 
ization j , 



Dj (q,a;„) = -- 



(q) 



(q)' 



(3.11) 



In Eqs. p.9p and p.lOp . q = p — p' and the symme- 



try relation \g°''^\ = \g \ for the coupling functions 
is used. Further, we neglect off-diagonal terms of the 
normal Green function. 

Equation p.lOp describes the normal contribution to 
the self-energy. The intravalley and intervalley scattering 
contributions enter on equal footing. In principle, the 
integrals in the right-hand side contain not only linear in 
e„ contributions, but also renormalize the Fermi energy 
and the spectrum. The latter types of the contributions 
do not lead to important changes of physical quantities 
and are neglected. 

Equation (|3.9p is the self-consistency equation for the 
order parameter Ap^^^. It is clear that the equation for 
the intravalley order parameter A™*™ — Ao,p,e„ (diag- 
onal elements of the matrices in Eq. p.9p ) differs from 
the one for the intervalley gap function A™*'^'" = Aq p 
(off diagonal elements of the matrices in Eq. (|3.9p ). For 
the former, only the intravalley interaction is important, 
while for the latter, both the intravalley and intervalley 
interactions contribute. It is clear that, provided the in- 
tervalley interaction is negative, the intervalley pairing is 
more favorable than the intravalley one. 

For explicit calculations, it is convenient to use the rep- 
resentation of the temperature Green functions in terms 
of retarded Green's functions 



dz 2ImF^(p, z) 
27r z — iSn 
dz 6j(q, z) 
2n z — ifln ' 



(3.12) 
(3.13) 



where is the retarded Green function, and bj = 
2lm.D^ is the phonon spectral function. 

Substituting Eqs. ((XT^ into Eq. we carry 

out the summation over the Matsubara frequencies and 
perform an analytic continuation iSn UJ + iS with an 
infinitesimal positive 5. 

Considering first the contribution of the intravalley 
pairing to the gap function (first line of Eq. p.9p ) we 
bring it to the form 



dp' 



dz r dz' 

TT- / TT- bj{p-p,z) 
27r . 1^ In 



. tanh^+coth^ ^^^,^^,^^ 

I i- I U! — Z — Z' + to 

dp' dz 



(2^)2 



2^ ypp^tanh— ImF«(p',z). 



(3.14) 



The contribution A^Jf" to the gap function coming from 
the intervalley pairing can be written similarly. 

Now, we will analyze the phonon and Coulomb parts 
separately. 



A. Phonon part 

In the case of large doping levels considered here one 
can use standard approximations well known in the the- 
ory of conventional metals. In particular, only momenta 
close to the Fermi surface can be taken into account. 
Reducing the dependence on the momenta p by the de- 
pendence on the unit vector n, p —ppn, we average the 
gap function Ap^g over the Fermi surface and introduce 
the quantity 



dn Ar 



(3.15) 
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where dn denotes the integral over all directions on 
the Fermi surface, and f is the density of states at the 
Fermi energy. The normalization is chosen in such a way 
that 



dn = 1. 



(3.16) 



The integral over the momentum in the right-hand side 
of Eq. (|3.14p reduces in this approximation to the form 



(...) 



V (fi) I dn 
Is 



d^pi-) (3.17) 



where v (p) is the density of states, Eq. (|2.13p . 

Using for the phonon Green function its bare value, 
Eq. (|3.1ip . such that 

6, (q,z) = 27r [(5(z + t^,(q)) - <5(z - u,iq))] (3.18) 

and integrating over the variable ^p, we write the phonon 
contribution I A™*''" ) to the gap function A™*''" of 
Eq. ((XTi)) as 

/CO 
Klliz.u) Im(F«(z))"dz. 
-OO 

(3.19) 

The intravalley phonon kernel K^\{z,lo) entering 
Eq. (|3.19p can be written as 



1 



Ctanh 
bj — z' — z -\- i8 



dz' ajjiz') 



I' coth 4p tanh — coth 4p \ 

X I 2Z 2T 2T |T 3 20 

UJ + z' - z + tS ^ ' 



where a^]^/(z) is the Eliashberg function for intravalley 
phonon scattering processes 

a?i/(2) - / dn dn'Y,\9lln,pHn-n'] 

J Sf ^ Sf j 

X (5(z-Wj(pF(n-n'))). (3.21) 



According to Eq. (|3.2ip this function contains the dou- 
ble average over the Fermi surface of the electron-phonon 
coupling function squared. 

The function F{z) in Eq. ([XTO)) equals 



F^d^, 



(3.22) 



where F^ is the retarded anomalous Green function ob- 
tained from the corresponding temperature Green func- 
tion, Eq. dXT]). 

This integration over ^ in Eq. (|3.22p results in a re- 
placement of the functions P± by P± given by 



P±(z) 



in sign(z) 



(3.23) 



Z(z)Jz2- |Ao,p,.±Aq,p,,| 



As concerns the off-diagonal parts of the gap function, 
we have to include the intervalley interaction processes 
into the self-consistency relation. Then we obtain 



(a: 



/ ph 



KUz,L.) Im(F«(z))''dz 



+Z-Huj) Kll{z,u)lu,{F\z)yUz, 

J —OO 

(3.24) 

with the intervalley kernel 
1 



m{z,u,) = -, 



dz' al^f{z') 



I tanh + coth ^ tanh — coth ^ 



z + i6 



uj + z' ~ z + i6 



(3.25) 



where the Eliashberg function for the intervalley scatter- 
ing equals 

"?2/(^) = / dn f^n'^ i |.9p^n,p^(n-n' 

J Sf J S f j i 

x5(z-a;j(pF(n-n')±Q)). (3.26) 

Similar calculations for the normal self-energy part lead 
to an expression applicable for the contribution of both 
inter- and intravalley phonon modes 

{l-Z{uj))uj (3.27) 
{Kll{z,u:)+Kll{z,u:))lTn{G''{z)) dz, 



where 



(3.28) 



The formulas derived in this section completely describe 
the effects of the electron-phonon interactions. Now, we 
will investigate the remaining parts of Eq. p.9p originat- 
ing from the Coulomb interaction. 



B. Coulomb part 

Calculating the Coulomb part in Eq. p.9p one should 
first renormalize the Coulomb interaction integrating out 
high energy degrees of freedo m^^'^"^ and thus reduce the 
original model to a model with a certain energy cutoff 
ujc, such that ujo uJc fJ-, where ujd = 2300 K. This 
renormalization is also logarithmic and can easily be car- 
ried out in the ladder approximation. The final results 
can be expressed in terms of the pseudopotentials U^^ 
and J7^2, respectively 



l + ^yiiln(iL)' 

i + :/(yii-hyi2)in(J^) 



(3.29) 
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where 



V 



a0 _ 



V 



Q/3 



PF(n-n') 



,sdn dn' 



(3.30) 



and the matrix elements V^^, V^'^ are given by (|2.151 

mni). 

Equations (|3.29p show that the effective Coulomb in- 
teraction C/q^ can not be very strong. Since the Coulomb 

potential decays as |q| ^ in the momentum space, the 
function V^^ is smaller than V^^ and the pseudopoten- 
tials do not differ much from each other. The pseudopo- 
tentials U^^ monotonically grow with increasing V"^ , 
but their values are limited by v^^. 

The renormalization of the Coulomb interaction p.29p 
is not important near the Dirac point because both v and 
fi are small, but can considerably reduce it in the region 
M ^ leV. 

As concerns the normal self-energy, a contribution 
coming from the Coulomb interaction is small and can 
be neglected)^ 

Using the pseudopotentials [/" and U^^ [Eqs. ([X^ ]. 
we write their contributions Ajl and A}? to the gap func- 
tion as 



[AL"*™)p= ^Z-^HC/" /tanh^Im(F«(z)) 



11 dz 
2^' 



(3.31) 

( A^O^ = (c) / tanh ^ Im {F^iz))'' ^. 



(3.32) 



Equations (|XTI?)) . ([X^ . (IX?T|) . and (j?:^ are basic 
equations of the electron-phonon theory of superconduc- 
tivity in graphene. Based on them, we can derive a for- 
mula for the critical temperature. 



C. Critical temperature 

In order to make an estimate for the critical temper- 
ature, we simplify our equations according to standard 
procedures. Our goal is to clarify which type of pairing 
is more favorable and to estimate the critical temperature 
rather than to calculate it from first principles. Within 
these procedures, the calculations become considerably 
simpler, but we believe that our goal is still achieved. 

Following Ref. [23|, we approximate the system of 
equations by linearizing their right-hand side with re- 
spect to the gap functions, and we approximate the 
phonon kernels Eqs. p.20p and (|3.25l) by the following 
expressions: 



tanh \z\ , \uj\ < ojd 



otherwise 



(3.33) 



where Xap, (a, /? = 1,2) are the intravalley and interval- 
ley electron-phonon coupling constants 



Aq/3 — 2 



«a/3/(^) 



(3.34) 



(Actually, we have to calculate only An and A12 because 
A22 = All and A12 — A21.) Equations (|3.33p show that 
the further calculations can be performed independently 
for the quantities with "11" and "12". 

The function Z{u}) coming from the normal self-energy, 
Eq. (|3.27p . is just a constant and can be written as 



Z=l + Aii + Ai2 = l + A. 



(3.35) 



Then, the fact that Z does not depend on frequency 
leads us to the conclusion that the gap function A^^ 
does not depend on the frequency for uj < ujjj either [see 
Eqs. (P3T|1 and 

By using Eqs. p.33p and linearizing the self- 
consistency equations (jXTO)) . ((X^ . (P3T|) . and ((X^ . 
we can derive equations for the critical temperatures of 
the intravalley and intervalley pairings J'^'"*''" and T"^^'^^, 
respectively. At the end, only the pairing with a higher 
critical temperature should be kept and used for the de- 
scription of the superconductivity. 

In order to obtain the equation for the critical temper- 
ature we choose the following form of the function A^: 



Aph UJ <ujd 

Ac <jJD < IjJ < (jJc 



(3.36) 



with constants Aph and Ac- 

Using the approximation, Eq. (j3.36p . we reduce the 
equation for the critical temperature J'™*''" to the form 



— tanh 

^ ^T- 



1 -f All A 



12 



All - Mil 



where the parameter fi^i equals 



Mil = 



l + i^Vii ln( 



(3.37) 



(3.38) 



Note that the approximation written in Eq. (I3.36P leads 
to the replacement of in the argument of the logarithm 
in Eq. by the Debye frequency w^) in Eq. ((X^ . 

The solution of Eq. (|3.37l) exists only when the right- 
hand side is positive. Therefore, the intravalley pairing 
is possible provided An > Mn- 

By calculating the integral over z in Eq. p.37p . we 
write the critical temperature T™*''" of the intravalley 
pairing explicitly 



1.13wd exp 



1 + A 
All - Mil 



(3.39) 



As concerns the intervalley pairing, both the intra- 
and intervalley phonon interactions contribute and we 
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come to the following equation for the critical tempera- 
ture T™*'^'' of the intervalley pairing 

^ tanh = y^n+>^i2 

with the renormalized Coulomb interaction given by 

HVii + V12) 



^12 



(3.41) 



The intervalley superconductivity is possible provided 
the condition A — An + A12 > /ii2 is fulfilled and we 
obtain the critical temperature T™*'^'' for this type of the 
superconductivity in the form 



rpl'. 



l.lSujo exp 



1 + A 



A 



12 



(3.42) 



As the constant /i^2 only slightly exceeds i^li for a 
strongly renormalized Coulomb interaction, the interval- 
ley pairing looks more favorable. Moreover, the electron- 
phonon coupling A entering Eq. p.42p is a quantity that 
can be extracted directly from the angle-resolved pho- 
toemission spectroscopy (ARPES), which simplifies esti- 
mates of the transition temperature T'™*'^'"^ Eq. (|3.42p . 
Clearly, the superconductivity is possible provided the 
condition 



A > fJ.*12 



(3.43) 



is fulfilled. Explicit estimates for the critical temperature 
rpmter ^^.g performed in Sec. V. We restrict ourselves with 
calculation of the temperature T™*'^'' because the critical 
temperature of the intravalley coupling T^"*^"-^ Eq. ()3.39p 
is lower than T™*'^'' for realistic parameters of and 
(1^2 ■ In addition, the intravalley pairing is sensitive to 
impurity scattering, which contrasts the intervalley pair- 
ing. The effect of the impurities on the two types of the 
superconducting pairings is considered in the next sec- 
tion. 



IV. IMPURITIES 

In the previous sections, we considered superconduc- 
tivity in clean systems. Usually, it is assumed that non- 
magnetic impurities do not affect the superconducting 
transition temperature. ""^^ However, the situation is not 
as simple for a system with several valleys, where some 
of the superconducting correlations can be sensitive to 
the impurities. We have considered the superconducting 
intervalley and intravalley pairing in the clean graphene, 
and now we will study effects of the potential impurities 
on these types of the superconductivity. 

In order to model this, we introduce an impurity 
Hamiltonian for the two-valley system. Generally, both 
intervalley and intravalley impurity scatterings are possi- 
ble. The most general form of the Hamiltonian for disor- 
dered graphene taking into account its Dirac-type spec- 
trum has been written (in absence of electron-electron 



interactions) in Ref. [26]. However, as we consider 
graphene for energies far away for the Dirac point, we 
introduce a standard impurity Hamiltonian in the mo- 
mentum space 



cr'-q.cr 



(4.1) 



p,q,o- 



with the momentum-dependent impurity potential u^. 
By using the isospin representation, we rewrite this ex- 
pression in the form 

2 

Himp = XI XI "q''«i,p+q,aa/3,p,<T, (4.2) 
a, 13 = 1 p,<i,cr 

where the functions u'^^ are related to the scattering po- 
tential as 



,11 



,22 



,12 



,21 



q+Q- 



Since intervalley scattering processes require a large 
momentum transfer, they can not be caused by Coulomb 
impurities of the substrate. On the other hand, vacancies 
in the graphene sheet, adatoms, surface ripples, or topo- 
logical defects can lead to both intravalley and intervalley 
scattering events.* 

For calculations, we use the standard diagrammatic ap- 
proach and treat the corrections in the Born approxima- 
tionii^i^l Studying the system far from the Dirac points, 
we consider only diagrams with non-crossing impurity 
lines. For the calculations, we assume that the effects of 
electroii-phonon and Coulomb interaction have already 
been taken into account according to Eqs. p.2l) . p.9l) . 
and (|3.10p . which determines Z„ and A. Calculating the 
corrections to these quantities arising from the impurity 
scattering we denote the renormalized values by Zn and 
A, respectively. 

By using the standard diagrammatic expansion in the 
approximation of non-crossing impurity lines, we obtain 
the modified self-energy equations 



dp' ,,n|2/^^p^.„ Fl?^^ 



(2^)2 I VV,e™ ^P 

dp' 



,12 I 



{Zn - Zn)ie 



(27r)2 I '1 I 



77121 

^P'.e„ 



nie. 





,121 







<^p',e™ 




(4.3) 



(4.4) 



where -Fp,£„ and Gp_£„ denote the renormalized Green 
functions. To obtain these functions, one has just to re- 
place Zn and A by Z„ and A in Eq. ((3^ . The further 
calculations are similar to those performed previously. 
We calculate the momentum integral in Eqs. (|4.3p and 
in the standard way and expand the right-hand sides 
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of the equations in the gap functions Ao,p, Aq^p, which 
gives us the possibihty to calculate the critical tempera- 
ture Tc- As before, the intervalley interactions affect only 
the intervalley gap and the normal self-energy. 
This leads to the following set of equations: 



1 



2t 



■^n^QjP — -^nAq p -I- 



4 + 


Aq,p + Ao,p 




1 




Aq,p 




2t / 








2 




Aq,p + 


Ao,p 




1 




Ao,p 





2rii 



(4.5) 



(4.6) 



(4.7) 



Aq^p + Ao,p 



Here, we have defined the elastic scattering time 



(4.8) 



In Eq. (|4.8D . th and T12 are intravalley and intervalley 
scattering times 



(4.9) 
(4.10) 



where riimp is the impurity concentration. Deriving 



'^11 — ^imp ^ 




dn 




T^2 — ^imp ^ 




dn 





Eqs. (|4.5I1 - (|4.10I) . we assumed as usual that the disorder 
is weak. Therefore, the main contribution in the integral 
over the momenta comes from the vicinity of the Fermi 
energy. 

Calculating T™*'='' we can put in Eqs. (|4.5M.7p Ao,p 
Ao,p = 0, which immediately leads to the relation 



A 



Q.p 



(4.11) 



because the normal and anomalous self energy renormal- 
izations of Aq.p cancel each other. Using Eq. (|4.1ip we 
conclude that both Aq^p and Aq^p must turn to zero at 
the same temperature T™*'^'' and this means the super- 
conducting transition temperature T'™*'^'' for the inter- 
valley pairing is not affected by the disorder. 

One can also come to this result by replacing the func- 
tions Zn and Aq p in Eqs. (j3.3p and (j3.7p by Z„ and 
Aq_p and using again Eqs. (|4.5p - (|4.7p . Then, one comes 
to Eqs. (|3.32p with Aq^p replaced by Aq^p, which leads 
to Eq. (EUl). 

On the other hand, we see that the cancelation of the 
normal and anomalous self-energies does not occur when 
calculating Ao,p at zero Aq_p, which indicates that im- 
purities influence this parameter. In fact, comparing this 
result with the conventional theory of paramagnetic im- 
purities in superconductorS)22i2^ we see that T12 plays the 
role of the scattering time on magnetic impurities. Thus, 
the intravalley superconductivity is completely destroyed 



3 
o 

g-0.10 

T3 
0) 
D. 
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FIG. 1. (Color online) Pseudopotentials /in (blue line) and 
Mi2 (green crosses) as function of the charge carrier density 
n. 



as soon as the inverse scattering intervalley time r-,^^ be- 
comes larger than the transition temperature T™*''° in 
the absence of the disorder. 

It follows from the results obtained in the present and 
previous sections that, by studying the possibility of su- 
perconductivity in graphene, it is sufficient to concentrate 
on the intervalley pairing. 



V. ESTIMATES 

Having derived the analytical expressions for the crit- 
ical temperature [Eqs. p.39p and (|3.42p ]. we should de- 
termine now the parameters A and fi* . Since we want to 
describe a graphene sheet where the Fermi level can be 
tuned, we examine the dependence of on the chemical 
potential fi. 

In order to estimate the Coulomb repulsion parameters 
^li and /z*2 [Eqs. (|3.38p and (|3.4ip ]. we use the screened 
Coulomb potential [Eq. (|2.12p ] and average this expres- 
sion over the Fermi surface in order to obtain vVn and 
i/Vi2 . These quantities allow us to calculate the Coulomb 
parameters fi^-^ and /i*2 as functions of the chemical po- 
tential /i using Eqs. p.38p and (|3.4ip . In order to be 
specific, we have chosen k — 2.5 (Ref. '23') for the value 
of the effective dielectric permeability of the substrate 
(occupying halfspace) entering Eq. (|2.10l) . This value 
corresponds to Si02 . The dependence of the parameters 
fill ^^"^ ^12 Oil the chemical potential fi (doping level) is 
represented for this value of k in Fig. [TJ Further, we use 
vq = 5.3 eVA for calculations. 

One can see from Fig. [T]that both pseudopotentials fj,l2 
and fill decay slightly with increasing electron density n. 

Calculation of the numerical values of the electron- 
phonon coupling constant A is more difficult because 
one has to know exact values of the matrix elements of 
the electron-phonon interaction. At the same time, the 
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electron-phonon coupling A determines the self-energy 
of the electron-phonon interaction and can be extracted 
from photoemission studies. Therefore, we simply take 
this value from literature. 

For a long time, there has been a rather poor agree- 
ment between theoretical results obtained using the local 
density approximation (LDA) and experimental values 
concerning the total coupling strength and the ratio be- 
tween the two nonequivalent coupling parameters. The 
electron-phonon coupling is also sensitive to the substrate 
(For details, see Ref. [29,] and citations therein). 

A possible source of the disagreement has been iden- 
tified in a recent paper,— where a copper substrate sub- 
stantially screening the electron-electron Coulomb inter- 
action was used and the agreement between the theo- 
retical results^ and the photoemission experiment was 
found. Both the theory and the experiment with the 
copper substrate lead to quite low values of the electron- 
phonon coupling constant A that remain below 0.05 for 
electron densities up to lO^^'cm^^. Using the metallic 
substrate one should assume that the dielectric perme- 
ability of the substrate k entering Eq. (|2.10p is a non- 
trivial function of the momentum. In order to avoid ad- 
ditional calculations for this system we note that even 
setting /i^2 = 0) the values A < 0.05 can not provide 
superconductivity with a noticeable transition tempera- 
ture. 

Measurements of the electron-phonon coupling in 
potassium-doped graphene on Ir(lll) substrateSS, have 
lead to the value A = 0.28 for a doping level of /i = 
1.29 eV (corresponding to the electron density n w 
1 X lO^'^cm"^). Such a value of A would lead to a rather 
high transition temperature. However, the authors of 
Ref. [s^l argue that the assumption of the linear spec- 
trum used in Ref. [29| leads to a considerably overesti- 
mation of A and expect lower values of this parameter 
corresponding to theoretical values of Ref. jsj . 

The authors of Ref. [sT'l suggest the following formula 
for the function A (yu) describing the dependence of the 
electron-phonon coupling on the chemical potential: 

A (/i) = 5.55CV^10~^cm (5.1) 

where n is the number of electrons per surface area de- 
pending on /i via /i = ^/^Tn and C — 1. 

This formula gives for n = 1 x lO^^cm^^ the value 
A = 0.056, which perfectly agrees with the experimental 
results of Ref. [30] for graphene on the metallic substrate. 

However, angle-resolved photoemission spectroscopy 
(ARPES) experiments^"— performed on doped 
graphene grown epitaxially on SiC lead to considerably 
higher values of A. Larger values of the coupling 
constants obtained for graphene on other substrates 
mean that the unscreened Coulomb interaction renor- 
malizes the electron-phonon interaction enhancing the 
latter. This conclusion correlates with the results of 
Ref. [23], where an enhancement of the intervalley 
electron-electron coupling constant A12 was predicted. 
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FIG. 2. (Color online) Critical temperature calculated with 
Eq. (|3.42p . The curves show Tc as function of the electron 
density n. Solid lines correspond to k = 2.5 while dashed lines 
are curves for k = 3.8. The bottom blue lines correspond to 
C — 3.5 and the top green lines to C = 5. 



So, we can try to use the values of the coupling constant 
A obtained for such a non-metallic substrate. 

According to a detailed analysis presented in Ref. [s^, 
the value of the coupling constant is 3.5 — 5 times larger 
than predicted theoretically^, which apparently implies 
that the coefRcient C in Eq. (|5.ip should take the values 
C ~ 3.5 — 5. The dielectric permeability of SiC equals 
K«3.8,^ 

As the constants C somewhat vary depending on the 
method of their calculation and the pseudopotential /i^2 
depends on the substrate, we simply draw in Fig. [5] the 
dependence of the critical temperature T"^*'^^ on the 
electron density n for several values of C and k using 
Eqs. ([CTI [5T|) . One can see from Fig. [H that 

the superconductivity is possible for realistic parameters 
characterizing the system and the transition tempera- 
ture T'*"*'^'' grows with increasing the electron density in 
graphene. Using the maximal possible value for C, Tc 
becomes very high reaching the value of 70 K for very 
high electron densities. This value of the critical tem- 
perature is apparently too high, otherwise it would have 
been observed in the experiment.^ Therefore, the value 
C — 5 does not look realistic. At the same time, the 
value C = 3.5 leads already to noticeable values of Tc. 

The analysis presented above was done assuming that 
the chemical potential ^ is far away from the Dirac points 
but is not close to the VHS. According to Ref. [32*], when 
approaching the VHS, the electron-phonon coupling A 
grows very fast, which would further increase the chances 
for the superconducting pairing. However, the linear- 
band estimation method used in Ref. [32] was shown to 
overstate the coupling,"^" and the growth of the coupling 
A near VHS obtained in the latter publication was very 
slow reaching the value A w 0.22. This value would still 
be sufficient for obtaining superconductivity with a rea- 
sonable critical temperature [see Eqs. p.42p ]. 
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Unfortunately, the experiment a^^'^'-'i'^^ — have not 
been supplemented by transport experiments on the same 
materials at low temperatures, and it is not clear whether 
the samples studied could be superconducting at low 
temperatures, or not. At the same time, the supercon- 
ductivity has not been observed^ in the transport mea- 
surements on graphene with Si02 substrate for the elec- 
tron density n up to 4x lO^'^cm"^ and more efforts should 
be expended to clarify the situation. 

VI. DISCUSSION 

In this work, we estimated the superconducting transi- 
tion temperature for graphene as a function of the chemi- 
cal potential fj, or area electron density n. We considered 
the case when the chemical potential is far away from the 
Dirac point, which corresponds to very high electron den- 
sity n. At the same time, we assumed that the chemical 
potential /i is not in the vicinity of the VHS. 

Starting with a model describing the electron-phonon 
and Coulomb interactions, we derived the Eliashberg 
equations for this system. Considering both anomalous 
and normal self-energies, we have obtained explicit for- 
mulas for the superconducting critical temperature that 
can be used not only for a weak electron-phonon cou- 
pling A, but also for A of order 1. We show that the 
Coulomb interaction in graphene is not very strong at 
high electron densities and does not necessarily destroy 
the superconducting pairing. 

As the parameters entering Eq. p.42p are not precisely 
known, we have drawn several curves in Fig. [5] corre- 
sponding to different values. It is clear that the critical 
temperature rather weakly depends on the dielectric per- 
meability of the substrate and other details characteriz- 
ing the Coulomb interaction. At the same time, the de- 
pendence on the electron-phonon coupling is strong, and 
we have shown curves corresponding to different values 
of these constants that may be considered as realistic. 

By estimating the pseudopotentials /i* describing the 
Coulomb interaction and using values of the electron- 
phonon coupling A extracted from photoemission exper- 
iments, we come to the conclusion that the transition 
temperature of the intervalley pairing J|J"*^'' can reach 
values exceeding 10 K for sufficiently high electron area 
density n. 

We have considered intra-and intervalley supercon- 
ducting pairing and demonstrated that the intervalley 
pairing is more favorable. Effect of disorder on the inter- 
valley superconductivity is weak but already a moderate 
concentration of impurities can destroy the intravalley 
pairing. All this means that the possibility of the in- 
travalley pairing can be discarded in realistic situations. 

According to previous findings, the coupling constant 
A can be considerably reduced provided a metallic sub- 



strate is used, which makes the superconductivity im- 
probable in such systems. 

Intercalating graphene by various materials may lead 
to an additional source of attraction between electrons 
and increase of the superconducting transition tempera- 
ture. 

The superconductivity becomes even more favorable 
when approaching the VHS. This is clear from the theo- 
retical point of view because the density of states diverges 
at this point, which should lead to a considerable increase 
of A. The region of electron densities of order 10^'^cm^^ is 
apparently already rather close to the VHS. This would 
imply that the approximation of the linear spectrum is 
no longer applicable. At the same time, the Fermi veloc- 
ity decreases when approaching the VHS, leading to an 
additional increase of the density of states and, hence, of 
the critical temperature Tc- 

A slow growth of the electron-phonon coupling near the 
VHS obtained in Ref. [35] using non-crossing self-energy 
diagrams indicates that this divergency is missed in this 
calculation. Moreover, using only non-crossing electron- 
phonon diagrams is not legitimate near the VHS and, 
therefore, the analysis of Ref. [35| is incomplete. 

Provided the electron-phonon interaction grows and 
the interaction remains essentially attractive one should 
expect at the VHS the conventional s-wave singlet su- 
perconductivity with a sufficiently high transition tem- 
perature. By analyzing logarithmically diverging dia- 
grams with the help of renormalization-group equations, 
a new type of unconventional (chiral) superconductivity 
was predicted recentl y^^'^^ in the situation when the in- 
teraction is repulsive. All this means that superconduc- 
tivity in graphene at high electron density is very proba- 
ble and we hope that it will be observed experimentally 
in the nearest future. 

Strictly speaking, superconductivity with the identi- 
cally zero resistance is not possible in 2D due to fluctua- 
tions of the order parameter and a flnite energy required 
for generation of vortices. The transition temperature 
Tc has been calculated in this work in the mean-field ap- 
proximation neglecting the fluctuations and vortices, and 
this is not justified. 

In practice, this means, however, that, instead of a 
sharp transition typical for 3D superconductors, one 
would observe a slower decay of the resistivity, which 
would make the transition rather broad. Although the 
resistivity does not become exactly zero in such a super- 
conducting state, its value can be extremely small and 
not distinguishable from zero in real experiments. 
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